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The observation of new physics events with large missing transverse energy at the LHC would 
potentially serve as evidence for the direct production dark matter. A crucial step toward verifying 
such evidence is the measurement of the would-be dark matter mass. If, for example, the invisible 
particles are found to have masses consistent with zero, it may prove very challenging to ascertain 
whether light dark matter or neutrinos are being observed. We assume that new invisible particles 
are pair-produced in a fi-like topology and use two jV/r2-based methods to measure the masses of 
the particles associated with the missing energy. We find that if the uncertainty associated with 
measuring the Mti endpoints is 0(1 GeV), the invisible particles must have masses greater than 
O(10 GeV) so they can be distinguished from massless ones. In general, the uncertainty associated 
with measuring the mass of the invisible particles increases as the mass decreases. If the recent 
results from the CoGeNT, DAMA/LIBRA, and CRESST experiments have indeed revealed the 
existence of light dark matter, our results suggest that it may be difficult for the LHC to distinguish 
dark matter from neutrinos solely via mass measurements. 

I. INTRODUCTION 

One of the main goals of high energy accelerator experiments like the LHC is the discovery, via direct production 
and subsequent detection of the associated decay products, of degrees of freedom beyond those present in the standard 
model. An advantage of direct production is that the various properties of the new particles, including their masses, 
can be measured, often with very good precision. Different methods exist for measuring the masses of the particles 
involved in different production and decay topologies [TJ [5] , but a unique problem arises when one or more of the final- 
state particles are "invisible," i.e., they do not interact directly with the detector. At a hadron collider, the existence 
of invisible particles can be inferred only via momentum conservation in the plane transverse to the beamline; only the 
vector sum of the transverse momenta associated with invisible particles can be reconstructed. If two missing particles 
in a single event each originate from identical decay chains, e.g., a fj-like topology 1 , the Mt2 variable [3] is a useful 
tool for extracting the masses of the parent, intermediate, and invisible particles in the decay chain. Remarkably, this 
includes the mass of the particles associated with the missing energy in the event [4j-[7|. 

Two popular M^ 2 -based methods for extracting the masses of the particles in a it-like decay topology are recon- 
structing the Mt2 kink j3HS] and the Mt2 subsystems method [TJ. The Mt2 kink method involves measuring Mti 
endpoints, which we call M™ 2 ax , for different values of an ansatz for the invisible particles' mass [JHB]. Analytical 
expressions can be fit to the distribution of M™ 2 ax as a function of the input ansatz, and the masses of the particles in 
the decay chain can be determined simultaneously from this fit. A kink exists in the M™ 2 ax distribution at the mass 
of the invisible particles, but as the mass of the invisible particles becomes light, it may be difficult to determine that 
the location of the kink is non-zero, given experimental uncertainties. Another Mj^-based method involves studying 
Mti subsystems and measuring three kinematic endpoints [7 . With the analytical expressions for the endpoints of 
each subsystem, one can simultaneously solve for the masses of the particles in the decay chain. If the mass of the 
invisible particle is light, the uncertainty for solving for its mass can become large. Because of the experimental 
uncertainties associated with these Mti methods, there is a minimum mass of the invisible particle above which it 
can distinguished from a massless particle at 95% C.L. 

Much of the literature to date uses SUSY models when addressing mass determination with Mti methods, taking 
the LSP to be 0(100 GeV) [5J El [5HUJ- The motivations for considering dark matter having these masses in the 
context of the MSSM and in light of experimental constraints are described in Ref. [18] and the references found 
therein. Recent data from the CoGeNT [TS], DAMA/LIBRA [HHH], and CRESST (22] experiments, on the other 
hand, hint at the existence of light dark matter (mass of order a few to 10 GeV) . Only a few examples exist in the 
literature which study light dark matter properties, including masses, at hadron colliders, e.g., (23U24| . In particular, 
if the dark matter is light, one may not experimentally be able to distinguish it from massless invisible particles, 
allowing models where neutrinos can fake dark matter signals at a hadron collider [25] . In our analysis, we estimate 



We define a tf-like event as any decay topology similar to the standard model (SM) leptonic decay of two top quarks: tt — > W + bW~b 
1+ v/bt~ V/h. We call the top quark the parent particle, the W boson the intermediate particle, and the neutrino is the invisible particle. 
All of these particles are considered to be potentially massive, while all other particles in the decay chain are assumed to be massless. 
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how massive the invisible particles must be such that they can be distinguished experimentally at a hadron collider 
from massless ones. 

The outline of our analysis is as follows: in Sectio n [IT we define the My 2 variable and express the previously- 
derived [3HS] expressions for its endpoints. In Sections III and IV we create pseudo-data and use the Mt2 kink and 
subsystems methods, respectively, to determine the mass that the invisible particle must have in order to distinguish 
it from a massless particle at 95% C.L. We do this for various masses of the parent and intermediate particles, and 
different uncertainties associated with the determination of the Mti endpoints. We conclude and discuss our results 
in Section IVl 



II. THE M T 2 VARIABLE 



Consider a general decay chain where a single massive particle, A, with mass rriA, decays to N massless, visible final- 
state particles and one potentially-massive, invisible particle with physical mass m. The individual three-momenta 
of the visible particles, Pi, are measured, and their four- momenta, p^ , are inferred with the massless approximation. 2 
The invariant mass of the visible system, M v { s , is defined as 

N 

MX = P^, P" = £>f. (1) 

i 

At a hadron collider, the sum of the transverse momentum of all final-state particles in the event is zero (to a good 
approximation), and the transverse momentum of the invisible particle, ft , can be inferred. We adopt an ansatz for 
the mass of the invisible particle, rh, and define the transverse mass variable Mt, the square of which is defined as 




A distribution of Mt values displays a "kinematic endpoint" or "edge" at the value of rriA when rh = m. 

Consider now an event where a pair of A's is created. Each parent particle and their daughters belong to a "branch" 
or "decay chain", where we add the label "(1)" and "(2)" to distinguish between the respective decay chains. The 
two parent particles, A^ 1 and A^ 2 \ eventually decay to and visible effectively-massless particles and a 

potentially-massive invisible particle with mass m. 3 The measured missing transverse momentum, ft , is the vector 
sum of the transverse momenta of the two invisible particles. Because there are two invisible particles in the final state, 
the Mt variable does not provide information concerning the masses of the particles in the decay chain. However, 
a generalized Mt variable, called Mti, can be introduced [3], and because we do not know how the transverse 
momentum is shared between the invisible particles, the square of Mt2 is defined as 



M T2 \ 



(rh) = min max \m t (p^ (1) , ; rh) , M T (p^ (2) , ; rh) 1 , (3) 



where and ' are free parameters over which the function is minimized, subject to the constraint kft = ft T - 

The distribution of Mt2 exhibits an endpoint, M™ 2 ax , at itia when rh = m. If rh ^ m, the Mti endpoint still exists, 
but its location does not equal rriA- In general, the change in the value of M™ 2 ax as a function of rh contains the 
information of the masses of the particles in the decay chains. A remarkable quality of the Mti variable is, for some 
topologies, the analytical form that describes M™ 2 ax when rh > m is different from that when rh < m. While these 
functions are continuous for all rh, their first derivatives are not, and a "kink" exists in the distribution of M™ 2 ax (m) 
at rh — m [H[6]. By analyzing the shape of the M™ 2 ax distribution as a function of rh, one can determine, in principle, 
the masses of the particles in the decay chains. 

For the rest of our discussions, we will assume a rt-like topology, as shown in Fig. [l] A^ — > B^D^\ followed 
by — ¥ C^E^\ where i = 1,2 and C is considered to be the potentially-massive invisible particle. We make 
no further assumptions regarding what particle types are D or E, only that they are effectively-massless and their 



2 The results presented here, of course, also apply for massive visible particles, as long as these are properly identified. 

3 It is unnecessary to provide any further details concerning the topology of the decay chain at this point. 
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FIG. 1: The decay topology of the pair production of parent particle, A, which decays into two visible effectively-massless 
final-state particles, D and E, and an invisible potentially-massive particle, C. Here, B is considered to be on-shell. 



three-momenta can be reconstructed. Particles A, B, and C have mass rriA, tub and mc, respectively, and we assign 
the ansatz fhc for the mass of the invisible particles. The analytical expressions for M™(mc) for the tt-like topology 
are, in the limit of no initial-state radiation (ISR) [51(7], 



M£f x (m c ) 



m A — m c 
2m 4 



2m A 



2m / 



fh c < mc, 



m A 
2 



2m j< 



m A 
2 



(4) 



m c , fhc > mc- 



By fitting the distribution of M™ 2 ax as a function of fhc with the expression for Af™ 2 ax (mc) in Eq. Q, one can 
simultaneously solve for mA, tub, and mc- 

Finding the kinematic endpoint of the Mt2 distribution can be difficult [26| . The endpoint feature can be obfuscated 
by decay widths, finite detector resolutions and lack of statistics around the kinematic endpoint. Additionally, it is 
often experimentally difficult to distinguish different types of final-state jets, e.g., if a reconstructed jet is due to a 
gluon or quark, so if there are jets in the final state of the pair-produced decay chain, the Mt2 endpoints can be 
contaminated with initial-state radiation jets [271128] . There can also exist combinatorial ambiguities associated with 
which jet is to be paired with which decay chain |29H31j . While these experimental issues are crucial and need to 
be ultimately understood in detail, we assume in this analysis that they only contribute to the overall uncertainty 
associated with determining the location of an Mt2 endpoint. This will become clear in the coming sections. 



III. M T 2 kink method 



The determination of M™ 2 ax (mc) for a given fhc will have an associated uncertainty, <je- To study the uncertainties 
associated with the determination of the mass of the invisible particle, one could consider producing Monte Carlo 
(MC) events of a tt-like topology with a detector simulation and using a procedure to measure the M™ 2 ax distribution 
as a function of fhc- Because of the significant amount of time it would take to create and analyze such MC samples, 
we choose to simulate simplified pseudo-data of the Mj^{mc) distribution. We take note that the individual values 
of M™2*{ m c) are highly correlated between values of fhc that are close together. This is easily understood since the 
events that populate the Mt2 kinematic endpoint for a given fhc are mostly the same, regardless of the choice of fhc- 4 
For simplicity, we assume that this correlation does not depend on the physical masses of the decay topology, 5 and we 
choose to sample M™(mc) in 0.25 GeV steps of fhc- We find that a smaller step size of fhc yields endpoints that 
are too correlated for adjacent values of fhc, and a larger step size would hinder the ability to resolve the position 
of the kink when it has a value close to zero. When simulating the simplified pseudo-data, one must take care to 
estimate the positive correlations between Mt2 endpoint measurements because ignoring to do so would imply infinite 
precision with which one could determine the M™ 2 a *(fhc) distribution and thus determine the masses of the particles 
in the decay chain with infinite resolution. 



4 The endpoint of the Mj-2 distribution is due to a particular momentum configuration of the final-state particles in the decay chains [5j E] . 

5 In general, highly-correlated values of f° r different values of mc involve a high percentage of the M™ 2 ax region from "unbalanced 
solutions" of Mt2 (the definition of which can be found in Refs |61 |32| V 
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The full summary of our study of the correlation between M.Ti{jhc) endpoints can be found in Appendix [A] We 
estimate a nearest-neighbor correlation (NNC) coefficient (as defined in Appendix [A| of 0.5 between neighboring 
M T 2{fhc) endpoints for 0.25 GeV steps in m C - 

With this choice of NNC, we generate fifty thousand pseudo-data M£? 2 ax (mc') distributions for different values of 
toc, fitting them with the analytical functions for Af™ 2 ax in Eq. Q, for which the physical masses of the decay topology 
are the fitting parameters. For each fitted distribution of M™ 2 ax , we histogram the fifty thousand best fit values of mc 
(which explicitly marginalizes over the uncertainties associated with the uia and mg fitting parameters), and from 
the width of this histogram, we estimate the 95% C.L. 6 associated with the value of mc- 

To draw similarities from a SM production of it events, i.e., for tua = 172 GeV, ms = 80.4 GeV (while still letting 
mc float), the uncertainties of the measurement of mc given a Ei = 1 GeV and a Ei = 5 GeV are shown in Fig. El 
From these figures, one can determine what the mass of mc must be in order to distinguish it from mc = at 95% 
C.L. If we call this value mg in , then mg in 9 GeV if a E = 1 GeV, and to™ 11 w 17 GeV if a E = 5 GeV. We also 
estimate m™ n when = 500 GeV and 77173 = 100 GeV, as shown in Fig.|3j and when 777^4 = 500 GeV and rag = 480 
GeV, as shown in Fig. [2] 




FIG. 2: The 95% C.L. for the mass of m c for m A = 172 GeV and m B = 80.4 GeV with (a) a Ez = 1 GeV and (b) a Ei = 5 
GeV, using the Mt2 kink method. The variable m™ ln is the value of mass of mc at which it can be distinguished from zero at 
95% C.L. 




FIG. 3: As in Fig. [2] for m A = 500 GeV and m B = 100 GeV. 



6 We use the term "C.L." to describe the uncertainty associated with the sampling method. 
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IV. M T 2 subsystems method 



To simultaneously measure wia, m B , an d met the M?2 subsystems method relies on the measurement of three 
independent endpoints: E 2W (M™ x (m c = 0) of the C 1 E 1 C 2 E 2 system), E 221 (M™ 2 ax (m c = 0) of the D 1 E 1 D 2 E 2 
system), and Ei m (endpoint of the invariant mass of the visible D 1 E 1 or D 2 E 2 systems). These endpoints can be 
expressed as functions of the physical masses in the decay chain. Following the naming convention in Ref. [7], 

E 221 = 1^4, (5) 
m A 



E 210 = - m & - TO ^ (6) 

V m 2 A 



/ (m\ - m|) (m| - m 2 c ) 

Eim ~ y ^| ■ (7) 

Given a measurement of the endpoints where i = 221, 210, or im, each with an associated uncertainty, a Bi , one 
can, in principle, invert Eqs. ([5|-([7| to solve for itia, mj, and mc as a function of the three E^s. However, because 
the measured endpoints may not have exactly the expected values, the system of equations may not be invertible 
without introducing large uncertainties. Instead of inverting the equations, we choose to perform a \ 2 fit using the 
measured values of Ei , , and the expected values of the endpoints as a function of the physical masses as expressed 
in Eqs. j5b-d7|). This x 2 function can be minimized with respect to m^, m Bl and m c , yielding the best estimates 
wiA, m %, an d m c for every set of Ei and <JEi- In the fit, we allow m 2 B and m 2 c to float negative but constrain itia to 
be positive. We assume that the uncertainties associated with the three E^s are uncorrelated. 

To estimate the uncertainties associated with our ability to measure mc, we generate fifty thousand pseudo- 
experiments, each with a set of three endpoint measurements, and each endpoint with the same uncertainty, (JEn 
for simplicity. A x 2 function is minimized for each pseudo-experiment, negative values of m c are set to zero (since 
they are unphysical), and the square root of the positive values of mc are histogramed. We find that the center of this 
distribution is centered about the physical mass mc and integrate about this center to find the 95% C.L. associated 
with the uncertainty of mc- 

In order to compare these results with those found using the Mt 2 kink method, we produce similar plots as in 
Figs. [2]|4j using the My 2 subsystems method. These are depicted in Figs. [5][7j We find very similar results for the 
value of m™ ln between the Mti subsystems and kink methods. This makes us confident that the procedure described 
in Appendix \K\ is sufficient to simulate the correlations between different Mt 2 endpoints as a function of mc- 

We estimate m™ m as a function of both and bib, as shown in Fig.[8j using the Mti subsystems method. We find 
when the masses of the parent and intermediate particles are close, the value of m™ n is higher as compared to when 
m>A and mj have a large mass difference. We choose not to compute the two-dimension plots of to™ 111 using the Mt 2 
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FIG. 5: The 95% C.L. for the mass of m c for m A = 172 GeV and m B = 80.4 GeV with (a) a El = 1 GeV and (b) a Bi = 5 
GeV, using the Mt2 subsystems method. The variable mJ5 ln is the value of mass of mc at which it can be distinguished from 
zero at 95% C.L. 





FIG. 7: As in Fig. [5] for m A = 500 GeV m B = 480 GeV. 
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kink method, as in Fig. [8j since creating the pseudo-data to make the M™ 2 ax (mc<) distributions is computationally 
expensive, and would require a sizable amount of running time. 




m A [GeV] m A [GeV] 

(a) (b) 



FIG. 8: The values of mJS ln , as a function of uia and tub, when (a) erg; = 1 GeV and (b) crgj = 5 GeV using the Mt2 
subsystems method. 



V. CONCLUSION 



The observation of new physics events with large missing transverse energy at the LHC would potentially serve as 
evidence for the production of dark matter. A crucial step toward verifying such evidence is the measurement of the 
masses of the would-be dark matter particles, i.e., the invisible particles. If, say, an excess of invisible particles is 
discovered and their masses are found to be consistent with zero, one can either conclude that there are new light 
massive invisible particles, their masses obfuscated by experimental resolution, or that there is new physics in the 
neutrino sector gjj]. As hinted by the CoGeNT Ql], DAM A/LIBRA [13 EI], and CRESST g3 experiments, dark 
matter may have a mass of O(10 GeV). Coincidently, assuming that dark matter is produced at hadron colliders as 
the decay product of a tt-like event topology and considering reasonable experimental uncertainties, we find that dark 
matter must have a mass greater than 0(10 GeV) such that it can be distinguished from neutrinos at 95% C.L. using 
Mr2-based methods. In general, the uncertainty associated with measuring the mass the invisible particles increases 
(decreases) as the mass decreases (increases). Our results suggest that, at the LHC, it will prove very challenging 
to distinguish light dark matter from neutrinos through mass measurements alone if the dark matter weighs around 
10 GeV or less. 

As seen in Fig. [8] the uncertainty associated with the measurement of the invisible particles' mass increases as the 
mass rng of the intermediate particles is close to the mass of the parent particles, and one finds the inverse effect 
if the intermediate particles are much lighter than the parent particles. In other words, a more precise measurement 
of the mass of the invisible particles can be made when the intermediate particles have a larger momentum in the 
center-of-mass frame of the two parent particles. Furthermore, the uncertainty associated with the measurement of 
the invisible particles' mass also increases as the parent particles' mass increases, for fixed mg/m/i. These features 
are due to the fact that the kink structure of the M™(mc) distribution becomes more pronounced when there is a 
large mass difference between the parent and intermediate particles. A similar scenario occurs when the intermediate 
particle is off-shell, such that mj> uia > m-Ci i- e -i the kink structure becomes more pronounced [7]. While a more 
pronounced kink structure will increase the precision at which one can measure the masses of particles involved in a 
given decay chain, it is not generally expected that it will change the shape of the uncertainty associated with mc as 
a function of the physical masses. 

The precision associated with measuring the mass of the invisible particles comes from how precisely one can 
locate the position of the Mj-2 endpoints with a model-independent procedure. Detector resolution, combinatorial 
ambiguities, and jet contamination are among the many significant obstacles for the precise determination of Mt2 
endpoints. We choose the la Gaussian uncertainties of <je — 1 GeV and <je = 5 GeV to be, respectively, a realistic 
estimate (compared to, say, reconstructed ti final states at the LHC) and a more conservative estimate for the 



uncertainties associated with the capabilities of the LHC [33 . If the uncertainty for locating an endpoint is greater 
than 5 GeV, we expect the precision with which one can measure the invisible particle mass to deteriorate accordingly. 

The individual events that populate the region about an M T2 endpoint for a given value of the ansatz mass are, to 
a good approximation, the same events that populate the endpoint for another nearby value of the ansatz mass. As 
shown in Fig. [9] (see Appendix |A|), the values of M™ 2 ax are not randomly distributed about the theoretical expectation, 
and while there are systematic uncertainties associated with Mt2 endpoint-finding procedure, there can be significant 
correlations between the measured values M™(mc) for adjacent values of rhc- For this reason, one cannot simply fit 
the measured value of Af™ 2 ax using Eq. Q while treating the individual values of M ™ 2 ax as independent measurements. 
In our analysis, we take this correlation into account when fitting with Eq. Q and marginalize over the other fitting 
parameters in order to correctly quote the uncertainty on the mass of the invisible particles. The majority of analyses 
in the literature ignore the correlation between values of M™ 2 ax for adjacent values of rh c , e.g., Refs. [51 151 [T21 ITTH2"4] . 
If these correlations are not taken into account when using the Mt2 kink method, we find that the uncertainties 
associated with the measured masses are underestimated. 

Other Mr2-based methods and topologies could have been considered in our analysis. All methods capable of 
extracting the mass of the invisible particles, however, take advantage of singularities in the kinematic phase space [341 
135] . Though some attention has been paid to other methods that can, in principle, provide more precise mass 
measurements [36j |37], it is likely that the limitations of the Mt2 kink and subsystems methods, discussed here 
in some detail, will be similar to the limitations of other methods as far as determining the mass of light invisible 
particles is concerned. It is also likely that the features associated with the change in the uncertainty of the invisible 
particle as the physical masses change are not a consequence of the M<r2-based methods we choose, but rather an 
inherent quality of measuring the mass of invisible particles. As already emphasized, the uncertainty associated with 
the mass of the invisible particle increases (decreases) as its mass decreases (increases), as can be seen in Figs. 2j[7 



In addition, as shown in Fig. [8j the uncertainty increases (decreases) as the difference between the mass of the parent 
and intermediate particles decrease (increase). We believe that these results are independent from the method used 
to extract the mass of the invisible particles. 
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Appendix A: Correlations between Mt2 endpoints 

To understand the correlation between endpoint measurements, we generate a MADGRAPH5 [38] sample of one 
hundred thousand SM pp — > ti events at y/s — 7 TeV, where the mass of the top quark is 172 GeV (corresponding 
to itia), and the mass of the W is 80.4 GeV (corresponding to ms). At the generator level, we fit the endpoint of a 
given Mt2 distribution with the following four-parameter probability density function: 

P (x;A,k L ,k R ,x ) = | Ae - kL(x - X( y a ™ x < x ^ (Al) 

where a; is a dummy variable for the x-axis of an Mt 2 distribution, A is a normalization constant, xq is taken to be 
the location of Af™ 2 ax , and k^ and kn are the exponential slopes on the left and right sides of the kink, respectively. 
We assume no combinatorial background associated with which jet is associated with which decay branch. The values 
of ^min and £ max are chosen a priori to be certain values about the expected value of M™ 2 ax . While picking values of 
x m i n and Xmax based off a theoretical expectation introduces a bias, we are only interested in extracting information 
concerning a rough estimation of how correlated adjacent endpoint measurements are as a function of rhc- For each 
Mt2{ttic) distribution, we perform an unbinned log-likelihood fit, where the log-likelihood function is defined using 
events with a value of Mf 2 between x m i n and 2: maX i 

ln£= ^ \np(xi;A,m L ,m R ,x ). (A2) 
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The extremization of ln£ is performed with MINUIT [33], and the results of the fits are shown in Fig.[9j The uncertainty 
associated with Xq is found by marginalizing over the uncertainties for A, mi,, and ijir and varying the log-likelihood 
about the minimum function by a value of 0.5 as a function of Xq. 
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FIG. 9: The fitted endpoints of Mti at the generator level with one hundred thousand events for tua = 172 GeV, raj = 80.4 
GeV, and mc = GeV using an unbinned endpoint fitting procedure. The dashed blue and dashed red curves are the theoretical 
expectations (the red dashed line is not shown since mc = GeV). Note that the points are not randomly distributed about 
the theoretical expectation. 

While this particular MADGRAPH5 sample vastly oversimplifies the type of distributions one would have to work with 
at a collider experiment, we only wish to extract information concerning the correlation between adjacent endpoints. 
To estimate the correlation between Mti endpoints and to create pseudo-data, we need to simulate a distribution 
that properly mimics the one in Fig. [9] while avoiding having to make thousands of MADGRAPH5 samples. To create 
pseudo-data that contain the correlations we see in Fig. [9j we need to quantify how much each endpoint measurement 
is correlated with others. 

We employ a simplified model of nearest-neighbor correlations. In this method, a single endpoint has its own 
uncertainty, CTj, and is positively correlated with the endpoint immediately to the left (except for the left-most 
endpoint, which is taken to have no correlations with any other endpoint). For simplicity, we consider that all 
measurements of M™ 2 ax have the same uncertainty, ue- We create a covariance matrix, which only includes nearest- 
neighbor correlations: 



Vu 



U E> 

NNC x cr|, 
0, 



* = J 

\i-j\=l 
otherwise, 



(A3) 



Here, NNC is the nearest-neighbor correlation factor. With this error matrix, one can iteratively generate pseudo-data 
that resemble a type of distribution we find with a MADGRAPH5 sample and a simple fitting procedure. We begin 
with the first point, no, at fhc = GeV, allowing it to be Gaussian-distributed about the theoretical expectation, 
/irj, with an error of <je- The next point, n 1; is Gaussian-distributed about the expectation with uncertainty oe 
but also is positively correlated to the the previous point. These endpoints are random variables sampled from the 
following probability density function (up to an overall normalization factor): 



,-(m-m) l<j h 



i = 



-(m-mf/al x e -(n i -fj,i)(m- 1 -fJ, i - 1 )V i /_ 1 j > 



(A4) 



Examples for the simplified pseudo-data for A/™ 2 ax are shown in Fig. 10 



While comparing the MADGRAPH5 and simplified pseudo-data plots by eye is suitable for our purposes, we wish 
to compare them more quantitatively. We can create a test statistic that estimates how much a given endpoint, rii, 
correlates with the endpoint immediately to its left, 



R 



(ni - Hi)(Tk-i - Mi-i) 



(A5) 



where ^ and are the theoretical expectations (given the physical masses of the system) for the M™ 2 ax distribution 
for bin i and the bin immediately to the left, If each endpoint is statistically uncorrelated with the one immediately 
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to the left, the distribution for R should be symmetric about zero. Fig. [Tl] shows a histogram of the values of R for 
the MADGRAPH5 sample and the probability density functions for pseudo-data for different values of NNC (these lines 
are determined by averaging over many generated pseudo-data). Guided by these distributions, we choose a value 
of NNC = 0.5 to perform the study discuss in Sec. |III| We choose this value because we interpret that the random 
jumps in Fig. [9]for < file < 10 GeV are innocuous qualities of our endpoint fitting procedure, and subsequently the 



histogramed values of —2 < R < in Fig. 11 for the MADGRAPH5 sample are not sampled from the true probability 
density function for R. Values larger than NNC = 0.5 where not studied because they lead to pseudo-data that 
deviated too much from the theoretical expectation, which is a feature not found in Fig. [9j Finally, we find that the 
correlation between adjacent values of rhc is not particularly sensitive to the physical masses of the particles in the 
decay chain, including ma- 



il) B. Gripaios, Int. J. Mod. Phys. A26, 4881 (2011), 1110.4502. 

[2] A. J. Barr and C. G. Lester, J. Phys. G37, 123001 (2010), 1004.2732. 

[3] C. Lester and D. Summers, Phys. Lett. B463, 99 (1999), hep-ph/9906349. 

[4] A. J. Barr, B. Gripaios, and C. G. Lester, JHEP 0802, 014 (2008), 0711.4008. 

[5] W. S. Cho, K. Choi, Y. G. Kim, and C. B. Park, Phys. Rev. Lett. 100, 171801 (2008), 0709.0288. 

[6] W. S. Cho, K. Choi, Y. G. Kim, and C B. Park, JHEP 0802, 035 (2008), 0711.4526. 

[7] M. Burns, K. Kong, K. T. Matchev, and M. Park, JHEP 0903, 143 (2009), 0810.5576. 

[8] C. B. Park, Phys. Rev. D84, 096001 (2011), 1106.6087. 

[9] K. Agashe, D. Kim, D. G. Walker, and L. Zhu, Phys. Rev. D84, 055020 (2011), 1012.4460. 

[10] M. Adeel Ajaib, T. Li, Q. Shafi, and K. Wang, JHEP 1101, 028 (2011), 1011.5518. 



11 



0.5r 
0.45! 

0.4! 
0.35! 

0.3! 
0.25! 

0.2! 
0.15! 

O.l! 
0.05! 

OS 



it 



MadGraph5 Simulation 
NNC = 
NNC = 0.4 
NNC = 0.5 



-1 



i-^v ■;■■!■■ 



1 2 3 4 5 
(n. - u.)(n M - \i j I ap lA 



FIG. 11: The variable R as defined in Eq. ( |A5[ ) for the MADGRAPH5 distribution as in Fig. [9] (black), and three values of the 
nearest-neighbor correlation factor, NNC. All distributions are normalized to unity. 
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